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ABSTRACT 

We have made the first map of CO(J=3-2) emission covering the disk of the edge-on 
galaxy, NGC 4631, which is known for its spectacular gaseous halo. The strongest 
emission, which we model with a Gaussian ring, occurs within a radius of 5 kpc. 
Weaker disk emission is detected out to radii of 12 kpc, the most extensive molecular 
component yet seen in this galaxy. From comparisons with infrared data, we find that 
CO(J=3-2) emission more closely follows the hot dust component, rather than the 
cold dust, consistent with it being a good tracer of star formation. The first maps of 
i?3_2/i-0i H2 mass surface density and SFE have been made for the inner 2.4 kpc 
radius region. Only 20% of the SF occurs in this region and excitation conditions are 
typical of galaxy disks, rather than of central starbursts. The SFE suggests long gas 
consumption timescales (> 10 9 yr). 

The velocity field is dominated by a steeply rising rotation curve in the region 
of the central molecular ring followed by a flatter curve in the disk. A very steep 
gradient in the rotation curve is observed at the nucleus, providing the first evidence 
for a central concentration of mass: Md yn — 5 x 10 7 M within a radius of 282 pc. The 
velocity field shows anomalous features indicating the presence of molecular outflows; 
one of them is associated with a previously observed CO(J=1-0) expanding shell. 
Consistent with these outflows is the presence of a thick (z up to 1.4 kpc) CO(J=3-2) 
disk. We suggest that the interaction between NGC 4631 and its companion(s) has 
agitated the disk and also initiated star formation which was likely higher in the past 
than it is now. These may be necessary conditions for seeing prominent halos. 

Key words: galaxies: individual (NGC 4631), galaxies: halos, ISM: bubbles, ISM: 
molecules, galaxies: ISM, ISM: structure 
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1 INTRODUCTION 

NGC 4631 (Fig. Q] Table EQ is an edge-orR galaxy that 
is known for a spectacular multi-phase halq_|. This galaxy 
is one of the targets of the James Clerk Maxwell Tele- 
scope (JCMT) Nearby Galaxies Lega cy Survey (NGLSjfl 
i|Wilson et al.ll2009l ; Iwturen et al.ll2010l ) whose goals include 
searching for molecular gas and dust in nearby galaxies and 
comparing the global properties of such systems. In addi- 
tion, the spatial and spectral coverage of the 325 - 375 GHz 
band presented by the new Heterodyne Array Receiver Pro- 
gramme - B-band (HARP-B) detector together with the 
wide-band Auto-Correlation Spectrometer Imaging System 
(ACSIS) has also made it possible to study individual galax- 
ies in the sample in some detail. For NGC 4631, our goals 
are to examine the CO(J=3-2) properties and distribution in 
this unique galaxy and to relate, where possible, the molec- 
ular emission to known outflow features. 

The strength, prevalence, and multi-phase aspects 
of the halo of NGC 4631 make this galaxy unique 
and an important target for disk-halo studies. The 
halo is observed in all ISM components, including cos- 
mic rays (CRs) and magnetic fields as indicated by 
polari z ed radio cont i nuum emission llGolla fc Hummell 
Il994al; iHummel et all Il99ll; IStrickland et all l2004al ). ~HI 



Randl Il994l; IRand fc Stonel Il99(ft 
Martin fc Kernl l200ll ; iRand et al.1 Il992l; lOtte et alj 



Hoopes et al.l 
Martin fc Ker 



dif fuse ionized gas 



19991) . dust dNeininger fc Dumkei 



2003 



1999 



200l|^ . molecular ga s (Rand 2000 a|), and hot, 

X-ray em itting gas fvogler fc Pietschlll996l;lTiillmann et all 



2006allbl; IWang"et all l200ll . 1 19951 : IStrickland et all l2004al : 
Yamasaki et al.ll2009lV Expanding shells h ave also been ob 



served in HI (|Rand fc van der Hulstlll993l ) and CO(J=1-0) 
|Randll2000al ). The halo extends over the entire star form- 
ing disk, reaching a variety of vertical heights, 2, depend- 
ing on th e component considered, frorrQ 900 pc in molec- 
ular gas jRandl l2000al) t o 10 kpc in the radio continuum 
i|Golla fc Humme]||l994ah . 

NGC 4631 is interacting with two companions, a dwarf 
elliptical, NGC 4627, 2.6' (6.8 kpc, AV ays = 64 km s^Jlto 
the north-west (Fig. [T]) and NGC 4656, another large edge- 
on galaxy about 32' (84 kpc, AV sy s = —40 km s _1 ) to the 
sout h-east, the result being 4 long intergalact i c HI s t ream - 
ers (|Weliachewl Il969l : I Weliachew et~aT1 1 19781 : iRandl Il994l ) 
stretching to ~ 42 kpc. The bases of these tidal streamers 
overlap with the halo of NGC 4631. In addit ion, three more 
faint companions have been detected in HI |Rand fc Stonel 
1 19961 : lRandlll994h as well as a faint optical dwarf galaxy 
candidate, NGC 4631 Dw A, 2.5 kpc below the plane of 
NGC 4631 (|Seth et al.1 l2005ah for which no redshift data 
are yet available. Of these companions, two (NGC 4627 and 
NGC 4631 Dw A) fall within the field shown in Fig. [1] 
and are marked with stars. Presumably, the star forma- 
tion (SF) activity (and hence the halo) in NGC 4631 has 



1 We take 'edge-on' to mean an inclination greater than 85° . 

2 The term, 'halo' is used to mean any extraplanar gas or dust, 
where we conservatively take 'extraplanar' to imply z ^ 1 kpc . 

3 http://www.jach.hawaii.edu/JCMT/surveys/ 

4 We adjust all values to a common adopted distance of 9 Mpc 
(Table 1), for comparative purposes. 



been triggered and/or enhanced by interactions. Interactions 
have also likely produced the observed thick stellar disk, 
i.e. the optical emission shows sech 2 (2/20) scale heights, 20, 
up to 1.4 kpc (depending on the stellar populat ion consid- 
ered) w ith detections to many scale heights in 2 (Set h et al.l 
l2005bl >. 

Studies of the mid-IR emissio n and dust proper - 
ties in NGC 4631 c a n be found i n iDraine et al.l d2007t) - 
iDumke et all (I2004D ISmith et al.1 d2007h , IStevens et al 



iBendo et al.1 i|2003l . 
l2007t ). Previous CO 



2006), and iDale et al 



observations in lower J 



transi tio ns have been carr ied out by Paglion e et al.l 
ll200ll), iGolla fc Wielebinskil ill994bl). [ Randl (I2000al1 . 
iTavlor fc Wand (|2003l ). and llsraell (|20od ). Ilsraell (|2009l ) 
also obtained a CO(J=3-2) measurement in a single beam 
at the center of the galaxy. Lim ited previous CO(J =3-2) 
mapping has been carried out bv lDumke et al.l (|200lT ) who 
detected emission only in the inner, 2.6' diameter region 
with a spatial resolution of ~ 23". The data presented here 
are of both higher resolution and sensitivity and, as will be 
shown, reveal the distribution of CO(J=3-2) both in the 
central regions as well as throughout the disk of NGC 4631. 

Sect. [2] outlines the observations and data reductions. 
In Sect. [3] we describe the CO(J=3-2) distribution, and 
will be considering the distribution in the disk and a com- 
parison with other wavebands, the CO(J=3-2) excitation, 
molecular mass and star formation, the velocity distribution 
and anomalous velocity and high latitude emission. Sect. [3] 
presents the discussion and Sect.[5]the conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Data were obtained of the 12 CO(J=3-2) (rest frequency, vq 
= 345.7959899 GHz) spectral line at the JCMT using the 
HARP-B front-end and the ACSIS back-end (Smith et al. 
2003). The HARP-B array contains 16 receivers in a 4 x 
4 pattern separated by 30" (about two beam- widths). In 
order to fully sample the field, the observing was set up as a 
sequence of raster scans in a "basket weave" pattern so that 
scanning was carried out in both the major and minor axis 
directions. The final sampling spacing was 1/2 of the beam 
si ze. Complete details of the observing set-up can be found 
in I Warren et al.l (|2010l ) and Table [5] In total, 14 scans were 
obtained over two nights under good conditions with the 
v = 225 GHz optical depth, T225GHZ, ranging from 0.051 
to 0.070 on January 5 and from 0.041 to 0.068 on January 6. 
The calibration sources were IRC+10216 and Mars. Pointing 
offsets over the course of the observations were 2" rms. 

Data reduction was initially carried out using the Joint 
Astronomy Centre version of the Starlink softwar^fl using 
the KAPPA, SMURF and CONVERT packages for edit- 
ing, cube making, and converting to flexible image transport 
system (FITS) format, respectively. Visualization programs 
such as GAIA and SPLAT allowed us to inspect the data. 
Two of the 16 receptors were not functioning properly and 
had to be removed from all data. The editing was itera- 
tive, beginning by inspecting cubes from each scan individ- 
ually, removing end channels, removing obvious interference 



See |http://starlink.jach.hawaii.edu| or lCurriel ((2008). 



CO(J=3-2) Observations of NGC 4631 3 




12 42 25 20 15 10 05 00 41 55 50 45 

RIGHT ASCENSION (J2000) 



Figure 1. CO(J=3-2) integrated intensity (zeroth moment) map of NGC 4631 superimposed on the Second Digitized Sky Survey (DSS2) 
blue image. The CO(J=3-2) field of view shown here is 9.8' X 4.1' in right ascension and declination, respectively. Contours are at 0.25, 
0.75, 2.0, 5.0, 8.0, 15, 24, and 45 K km s — 1 and the peak value is 57.7 K km s —1 . The moment-generating routine applied a spatial 
gaussian smoothing function of 40" FWHM and a cutoff level of 0.018 K before integrating over the unsmoothed cube from 437 km 
s~ 1 to 821 km s — 1 (see text). Red stars denote the centers of NGC 4631 (Table [TJ, the dwarf elliptical companion, NGC 4627 to the 
north-west, and the optical dwarf galaxy candidate, NGC 4631 Dw A, to the south-west. The red plus (+) and two crosses (X) denote 
the approximate centers of the CO(J=1-0) shell found by Rand (2000a) and the two HI supcrshells found by Rand & van der Hulst 
(1993, their Fig. 3), respectively. Features discussed in Sect. QTTJare labelled. 



spikes, binning to 20 km s _1 , removing a linear baseline, 
and then collapsing the cube to inspect the total intensity 
(zeroth moment) map. This usually revealed pixels that had 
obviously poor baselines for any given scan. Poor data points 
were removed from the unbinned, unbaselined data, and the 
process repeated, as required. All scans of the edited, but 
otherwise original resolution data were then combined into 
a single cube, using a 'SincSinc' kerneQ with a cell size of 
3.638 arcsec (1/4 of the beam) and then saved in FITS for- 
mat. 

The fits cube was then read into the Astronomical Im- 
age Processing System (AIPS) package for the remainder of 
the processing and analysis. The data were box-car binned 
to a velocity resolution of 10.4 km s~ and a linear base- 
line was removed, fitted pixel by pixel. Some further minor 
editing was then carried out in AIPS. In addition, residual 
baseline curvature was still evident in some sections of the 
cube and these were then flattened using a 3rd order poly- 
nomial. 

The final data were subsequently corrected for the main 
beam efficiency, t)mb = 0.60 (estimated uncertainty between 



7 The MAKECUBE routine in SMURF was used. 



10 and 15%) in order to convert into units of main beam 
brightness temperature, Tmb- Final channel maps for those 
channels that display emission are shown in Fig. [2^, and b. 
The resulting measured rms noise (Table[5J) per channel, met 
our goal for the NGLS^. An examination of each individual 
HARP beam from independent observations of Mars shows 
that sidelobes of order 3% of the peak occur at distances of 
24" from the beam center which is well within our estimated 
uncertainties (see below). Sidelobes at larger distances are 
estimated to be less than this (P. Friberg, private commu- 
nication) and therefore negligible. 

A total intensity (zeroth moment) map was then made, 
which involved smoothing the data cube spatially, imposing 
a flux cutoff based on the smoothed cube, and then integrat- 
ing in velocity over the original cube, using only those pixels 
that, in the smoothed cube, were above the flux cutoff. The 
result is shown in Fig. [TJ with cutoff and smoothing details 
given in the caption. The field of view has also been reduced 



8 The target rms per 20 km s _1 channel was 0.030 K (Tmb) 
which matches our measured value with binning to the wider 
channel. Note, however, that the noise increases towards the map 
edges. 
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Table 1. Basic Parameters of NGC 4631 



Table 2. Observing & Map Parameters of NGC 4631 



parameter 



value 



parameter 



value 



Hubble type_^ 
RA (J2000)~h m s) b 
DEC (J2000) (° ' "£ 
Vsys (km s- 1 ^, 
D (Mpc]£ 
a x b (' x 
PAf 
i9_~ 

Lfir (erg s 
Lfir/DI 
SFK FIR (Mq yr" 1 ^ 
SFR Ha , corr (M Q yr- 1 ^ 

(M £ 
M dust (M ^ 



i/Dfg (10 40 erg s" 



kpc" 2 ^ 



SB(s)d sp or Sc sp 

12 42 08.01 

32 32 29.4 

606 

9.0 

15.5 x 2.7 



9.5 x 10 43 

7.76 

4.3 

2.4 

1.0 x 10 10 
9.7 x 10 7 



a iButa et al.l J2007t) . 

" IR center at 2^tm from the Nasa Extragalactic Database (NED). 
c Systemic HI velocity (heliocentric, optical definition) from NED. 
d Distance (e.g. iKennicutt et aljf2003h . 

e Optical major X minor axis, me asured to the 25 blue mag 
per square arcsec brightness level llde Vaucouleurs et al.|[l99ll) . 
f Position angle |de Vaucoule urs et aTlll99l[) . 
9 Adopted inclination (e.g. Ilsraelll2009l y 

Far-IR luminosity flTiillmann et al.l l2006a). adjusted to 
D = 9. Mpc. 

' From Friillmann et al ] l|2006ah . where D25 is the galaxy's 

diameter measured at the 25 blue mag per square arcsec level. 

3 Far-IR SFR iTiillmann et al.ll2006a| y See also Sect. [331 

^ Ha SFR corrected for extinction according to the formula of 

Calzetti et al. [j2007t) , but using the Ha correction of 

Zhu et all <200S|) ; t he Ha luminosity of [Hoopes et alj <199St) 

and A 24 fim flux of lDale et al. I fcOOTf ) have been used. 



HI mass (Rand 1994) 



Dust mass (Bcnd o et alj|2006h . 



slightly (to 80% of the original) in order to trim noisy edge 
points. 

The noise, taken to represent the uncertainty due to 
random errors in our total intensity map, is ~ 1.6 K km 
s _1 , based on the noise per channel and number of chan- 
nels entering into the sum at any typical position contain- 
ing emission. In addition, we can compare our map to one 
obtained independently from the same data but using differ- 
ent software, e diting, pixel size, ch annel width, and baseline 
smoothing (see lWarren et all feoiO). The maxima of the two 
maps differ by 4%, and a histogram of the difference map is 
well fit by a Gaussian with a peak at 0.20 K km s _ and a 
standard deviation of 1.0 K km s" 1 . These differences, the 
largest of which are likely due to differences in baseline flat- 
tening, are less than our estimated uncertainty above. Aside 
from these random errors, an absolute calibration error of 10 
to 15% is present due to uncertainties in the value of 7]mb', 
this uncertainty affects all pixels and does not change the 
appearance of the map. We have also compared the inte- 
grated i ntensity from the center of Fig. [T] to the single value 
given bv llsraeli (|2009h who also used the JCMT but with a 



Observing Date Jan. 05 & 06, 2008 

Total bandwidth 1 GHz 

Original channel width 0.488 MHz (0.43 kms" 1 ) 

Velocity-binned channel width 11.7 MHz (10.4 km s — 1 ) 

Angular resolution 11 14.5" 



rms (Tmb) P er channel'' 



0.034 K 



a Average of HARP beams. 

" For 10.4 km s -1 channel width. 

14.5" beam agrees with his result of 22 ± 3 K km s _1 in a 
14" beam at the same central position. 

The ancillary data used in this paper were taken from 
the Spitzer Infrared Nearby Galaxies Survey (SINGS) Ancil- 
lary Data Archival unless otherwise indicated. The archive 
contains primarily Spitzer data but also includes ancillary 
images such as the Ha image used in this paper. More 
information about the SINGS program can be found in 
IKennicutt et al l (|2003l ). 



3 RESULTS 

3.1 The CO(J=3-2) Distribution 

As shown in the channel maps (Fig. [2]) the east side of the 
galaxy is receding and the west side is approaching. Both 
the channel maps and the total intensity map (Fig. [TJ show 
that the strongest emission is concentrated in a region of 
diameter, 3.8' (10.0 kpc) centered on the nucleus (labelled 
the 'central molecular ring' in Fig. [T] for reasons outlined in 
the next subsection). At larger radii, there is weaker more 
extended disk emission. Several other features which we re- 
fer to below are also labelled in Fig. \T\ Fig. shows a slice 
in emission along the major axis of Fig. [1] at a position angle 
(83.5°) chosen so that it passes through the two broad max- 
ima on either side of the nucleus (rather than the global op- 
tical major axis position angle of Table [1]). Structures along 
the CO(J=3-2) major axis are well represented by this slice. 
In the next subsections, we discuss features associated with 
the disk of NGC 4631; discussion of the vertical distribution 
is deferred until Sect. 13.71 



different receiver. Our mean value of 24 ± 3 K km s in a 



3.1.1 The Strong Central Molecular Emission 

The strong central molecular emission is characterized by 
two peaks on either side of the nucleus with a central mini- 
mum between them, and a slight major axis curvature that 
is well-known in this galaxy. This central emission extends 
between a minimum in the emission on the east (the 'east- 
ern minimum') and a gap in the emission on the west (the 
'western gap'), as labelled in Fig.[T] The eastern peak is 15% 
higher than the western on e, in agreeme n t wit h the rudi- 
mentary CO(J=3-2) ma p of Dumke et al. If200j). Since the 
CO(J=1-0) distribution (|Golla fc Wielebinskj|l994b1 : iRandl 



http://irsa.ipac.caltech.edu/data/SPITZER/SINGS 
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Figure 2. (a) CO(J=3-2) channel maps. Contours are at 0.10 (3<r), 0.20, 0.30, 0.50, and 0.75 K. The declination and right ascension 
scales are shown in Fig.[2p. The velocity of each channel is given at upper right and the beam is shown at the lower left of the first frame. 



2000a) also shows a central minimum, the observed central 
CO(J=3-2) minimum cannot be a result of lack of sufficient 
excitation but must reflect a true minimum in the molecular 
gas distribution between the two peaks. A similar structure 
is also observed in dust emission (e. g. iBendo et all 1 2002; 



iDumke et all |2004| ; IBendo et all 120061 ; IStevens et al 



2005); 



we discuss comparisons with other wavebands in Sect. 13.21 
Fig. |3] clearly shows that this double-peaked central 
molecular emission dominates the CO(J=3-2) distribution. 
The red curve represents an edge-on Gaussian ring with peak 



amplitudes that are slightly asymmetric and whose fitted pa- 
rameters are given in Ta ble [3] The em ission was modelled 
follow ing the method of llrwinl (|l994h and llrwin fe Sofud 
( 1996) which reproduces the line profiles of the cube. We 
imposed a cutoff radius at 100" in order to model the strong 
central emission only. The displayed model profile was then 



obtained in the same way as the data slica I. This Gaus- 



The model is symmetric but the amplitude was arbitrarily 
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Figure 2. (b) CO(J=3-2) channel maps as in Fig. [2^. 
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Table 3. Parameters of Modelled Central Molecular Ring 



parameter 



value 



RA (J2000) (h m 

DEC (J2000) (° ' ") a 

i (deg£ 

Ro (", kpc£ 

Do (", kpci£ 

Di (", kpc£ 



12 42 7.7 ± 0.4 
32 32 30 ± 5 
89 ± 4 

42 ± 3, 1.8 ± 0.1 
6.4 ± 0.7, 0.28 ± 0.03 
2.1 ± 0.7, 0.09 ± 0.03 



a Ring center position. Uncertainties indicate the variation 
that produces an estimated increase of la to the residuals. 



Best fit inclination. 



c Galactocentric radius of ring peak. 
^ Outer Gaussian scale length, i.e. n(r) 



i.e. n(r) = no exp(—r 2 / (2D%)) 
where n(r) is an in-plane density, no is the density at Rq, and r 
is a radial distance measured outwards from Ro. 
e Inner Gaussian scale length, as in Footnote [5] with Di 
replacing D and r measured radially inwards from R . 




100 -100 

Major Axis Distance (arcsec) 



sian ring describes the central emission very well except for 
excess emission in the wings at larger radii and also some 
departures near the nucleus. The inset shows the residuals. 

The geometry of the central molecular gas distribu- 
tion may be more complex than a simple ring. For exam- 
ple, bright spiral arms, rich in molecular g as of the kind 
seen in more face-on galaxies such as M 51 l)Brunner et al.l 
2008), may mimic a ring or pseudo-ring when observed edge- 
on. There has also been some suggestion that NGC 4631 is 
barred. For example, there has been uncertainty in th e op- 
tical classification (see Table [TJ and iRov et al.l l|l99ll) have 
suggested that the entire region that we have modelled as a 
ring could be a large bar. A molecular bar could also mimic 
the distribution shown in Fig. [3] provided its density distri- 
bution is peaked at the bar ends. W hile this is possible, it has 
been shown (see iKuno et al.l 120071 . and references therein) 
that barred galaxies are much more likely to show strong 
central peaks or concentrations of CO, rather than t he cen - 
tral minimum that we observe in NGC 4631. iRandl (|2000al ) 
also finds no evidence for a bar, and a new high quality opti- 
cal mosaic in gri colours from the Sloan Digital Sky Survey 
shows no clear bar in this galaxy eithej"1. 

In summary, other specific geometries could be invoked 
to explain the strong central molecular emission, but since 
the Gaussian ring can do so with very few free parame- 
ters, we will use the term 'central molecular ring' to de- 
scribe this region. Its velocity distribution will be discussed 
in Sect. 13.6.11 but it is worth noting here that the central 
molecular emission is kinematically distinct from the outer 
emission. 



Figure 3. Slice along the major axis of the total intensity 
CO(J=3-2) distribution shown in Fig. [T] averaged over a width 
11" (3 pixels) in z. East is on the left and west is on the right. 
Offsets are with respect to the infrared center (Table [TJ. The thin 
red curve shows a gaussian ring with parameters as given in Ta- 
blc[3] The inset shows the residuals (data - model) for the central 
molecular ring. 

3.1.2 The Nucleus 

As pointed out by iGolla fc Wielebinskil l|l994bl ) , there are 
disagreements in the location of the center of the galaxy, 
depending on how the center is defined or which compo- 
nent is considered, consistent with a system that is highly 
disturbed. Fig. [TJ for example, shows a central minimum in 
the CO(J=3-2) distribution (marked with a 'C'). This mini- 
mum agrees with the minimum seen in th e CO(J=1-0) map 
l|Golla fc Wielebinski|[l994bl ; lRandll2000al ) but is 15" to the 
west of the infrared (IR) center of Table [TJ (marked with a 
star). 

The center of our modelled ring (Table [3J, however, 
agrees with the IR center within uncertainties, but not with 
the positio n of C. The IR center coincides with the central 
radio peak l|Gollalll999l ). and our global CO(J=3-2) flux (not 
just the ring) is more symmetric about the infrared center 
than about C. In Sect. 13.61 we will also present a dynamical 
argument for the nucleus to be more closely represented by 
the IR center. Therefore in this paper, when we refer to the 
center or nucleus, we mean the IR center of Table QJ Note 
that there is a small peak in CO(J=3-2) right at the nucleus 
as can be seen in Fig. [3] (see also Sect. I3.6TTI) . 



fitted and the slight asymmetry was reproduced with an insignif- 
icant adjustment of the position angle. The results of Table[3]are 
not affected by this. 

11 See http://cosmo.nyu.edu/hogg/rc3/ courtesy of David W. 
Hogg, Michael R. Blanton, and the Sloan Digital Sky Survey Col- 
laboration. 



3.1.3 The Weaker, Extended Disk Emission 

In addition to the strong central molecular ring, we also 
observe weaker CO(J=3-2) emission at much larger radii, 
evidently associated with the larger rotating galactic disk 
(for example, see velocities of 499.3 km s _1 and 738.3 km 
s _1 in Fig. [2}. This weaker broader disk emission appears 
distinct from the central molecular ring, separated from it 
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by the eastern minimum and western gap (see Fig. [TJ; at 
these two locations, there is also an abrupt change in the 
rotation curve gradient (see Sect. 13. 6|) . We will refer to this 
emission as the 'outer disk' to distinguish it from the central 
molecular ring. 

This outer CO(J=3-2) disk tends to follow the optical 
disk emission. For example, the far eastern emission, cen- 
tered at a right ascension of about 12 h 42 m 21 s (the 'b ump , 
Fig. [1} maintains the bulging shape of the underlying op- 
tical disk. On the far western side there is a distinct cor- 
rugation in the CO(J=3-2) emission centered at a right as- 
cension of « 12 h 41 m 52 s (the 'kink', Fig. Q}. The outer 
dis k also harbours the two expan ding HI supershells found 
by iRand fc van der HulstJ l| 19931 ) (marked with crosses m 
Fig. D. 

The maximum extent of the detected emission is 3.5' 
(9.25 kpc) to the east of the nucleus and 4.7' (12.4 kpc) to 
the west. This radial extent exceeds that of previous CO ob- 
servations in any transition, although t he CO (J=1-0) emis- 
sion detected by iGolla fc Wielebinskil ljl994bh is almost as 
extensive to the west. 



3.2 Comparison of Disk Emission with other 
Wavebands 

The CO(J=3-2) emission is displayed in a number of over- 
lays in Fig. [4] which show distributions that are taken to be 
proxies for unobscured star formation (Ha), the hot dust 
distribution (MIPS A 24 /im) and the cold dust distribution 
(MIPS A 160 /Ltm). For the high resolution, high contrast A 24 
/im map, we also show an inset of the central region. It is 
well known that molecular gas, dust and star formation cor- 
relate and we can see evidence for this in these overlays. For 
example, the CO(J=3-2) bump on the eastern side is evident 
in the Ha map as well as the A 24 /im map, and the gen- 
eral trend of strong central emission and weaker secondary 
peaks at larger radii is common to the CO(J=3-2) and the 
two dust-emitting bands. 

To explore these relations more quantitatively, Fig. [5] 
further shows comparative major axis slices of CO(J=3-2), 
A 24 /im, A 160 /im, and Ha emission at a common spatial 
resolutiorQ and normalized to their peak flux values. Note 
that these slices contain most of the emission in any of the 
maps, since the width of the slice is 40 arcsec. Here, the gen- 
eral trend of the two dust components following the molec- 
ular gas distribution is again seen, but the departure of the 
Ha emission is more obvious since this latter component is 
strongly affected by dust obscuration. The correlation co- 
efficient is 0.99 between CO(J=3-2) and A 24 /im emission 
as well as between CO(J=3-2) and A 160 /im, whereas it is 
-0.001 between CO(J=3-2) and Ha emission. 

Focussing only on molecular gas and dust, Fig. [5] also 
shows that the CO(J=3-2) emission follows the hot dust 
emission (A 24 /im) more closely than the cold dust (A 160 
/im). Colder dust displays a broader distribution that de- 
clines more slowly from the central molecular ring. Between 



For the Spitzer images, resolutions were matched using con- 
volution kernels applicable to both th e beam size and shape (see 
iGordon et al"1l200Sl ; iBendo et ai1l2010l) . 
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Figure 5. Comparison of major axis slices at 40" resolution, av- 
eraged over 40" in z and normalized to their peak values. The 
positive x axis is to the east. The thick red curve shows the 
CO(J=3-2) emission taken from Fig. [l] and the thin red flanking 
curves show the extent of its error bars (which are the dominant 
uncertainties in the figure). The steep fall-off around +200 and 
-250 arcsec is due to the moment-generating routine. The black 
curve is the A 24 /im emission from Fig. [4b, the medium grey 
curve represents the A 160 /im emission from Fig. [4ji and the light 
grey curve represents Ha from Fig. [4^. 

± 200 arcsec, for example, the mean of the ratios, CO(J=3- 
2)/A24 /im and CO(J=3-2)/A 160 /im along the major axis 
(not shown) are 0.91 and 0.35, respectively, supporting the 
closer connection between CO(J=3-2) and hot dust. This 
result is consistent with the fact (see next section) that 
CO(J=3-2) is a good tracer of star formation and we would 
therefore expect it to be more closely aligned with hotter 
dust. Cold dust is more likely to trace both molecular gas 
farther from star forming regions as well as the more broadly 
dist ributed atom i c com ponent . 

IBendo et al.l l|2006l ) have found that there is no signifi- 
cant variation with radius between the A 70 /im, A 160 /im, 
and A 450 /im emission in this galaxy and so we would ex- 
pect plots of A 70 /xm and A 450 /im emission to follow the 
cold dust emission of A 160 /im displayed in Fig. [5] Available 
A 450 /im data do not have sufficient signal-to-noise to test 
this but the close relation between A 70 (im and A 160 (im 
data allows us to repeat the above analysis at a higher spa- 
tial resolution (18") using A 70 /im as a cold dust indicator. 
Again we find the same conclusion that CO(J=3-2) more 
closely follows hot rather than cold dust. 

Fig. [S] appears to show an exception to this result right 
at the nucleus where there is a peak in the A 24 /im hot dust 
distribution but minima in both CO(J=3-2) and A 160 /im 
cold dust. However, at higher spatial resolution (see the inset 
of Fig. |4p) we see that the nuclear peak at A 24 /im is due to 
a strong 'hot spot' located within a region of approximately 
17 arcsec (740 pc) diameter. As pointed out in Sect. 13.1.21 
there is also a small CO(J=3-2) peak at the nucleus so the 
CO(J=3-2)/hot dust relation appears to hold even there, 



CO(J=3-2) Observations of NGC 4631 9 




12 42 25 20 15 10 05 00 41 55 50 45 




12 42 25 20 15 10 05 00 41 55 50 45 

RIGHT ASCENSION (J2000) 



Figure 4. Overlays of the lowest and second highest contours of CO(J=3-2) emission from Fig.[T]on greyscales of (a) the Ha image, (b) 
the MIPS 24 fim image, and (c) the MIPS 160 (im image. The greyscale images are shown at their original resolutions whic h are 3.0" for 
the H a image (from measurements of point sources in the field), 6" for the 24 /im image, and 40" for the 160 ^m image jBendo et al.l 
2006). The greyscale ranges are arbitrary and have been chosen to emphasize low to moderate intensity features. For the high contrast, 
high resolution 24 fim image, we also show the central region in an inset. The star denotes the infrared center. 



although relative emission strengths may vary. We will show 
in the next section that the star formation rate is higher at 
the nucleus than in the immediately surrounding region. 



3.3 CO(J=3-2) Excitation 

We can study the CO(J=3-2) excitation in NGC 4631 by 
forming a ratio map of CO(J=3-2) to CO(J=1-0) emission. 
To this end, we use CO(J=1-0) data obtained using the 
Berkeley Illinois Maryland Array (BIM A) , ori ginally at reso - 
lution of w 8", kindly supplied by R. J. Rand (iRandlbOOOaT ). 
The total intensity CO(1-0) map from lRandl |2000al ) (his fig- 
ure 1) was interpolated to the same grid as our total intensity 
CO(J=3-2) map and then both images were smoothed to 17" 
resolution to reduce noise in the ratio map. The CO(J=1-0) 



map was then converte to units of K km s 1 . As a check 
on the CO(J=1-0) map, we made comparisons with previ- 
ous singl e dish CO(J= 1-0) data from the literature. As indi- 
ca ted bvlRand ll2000ah. the BIMA total flux agrees with that 
of lGolla fc Wielebinskil (|l994bh within uncertainties. Also, 
the BIMA integrated intensity agrees with the Institut de 
Radioastronomie Millimetrique (IRAM) va lue at the sa me 
central position and resolution as given in Israel (20090 
A comparison of major axis slices between the BIMA data 



13 The Planck relation was used. In the Rayleigh-Jeans limit, the 
conversion, for a beam solid angle of f2f, = 7.7 X 10~ 9 sr, is 1 Jy 
beam -1 km s -1 = 3.1 K km s . 

14 The BIMA result is 41 K km s _1 compared to the the IRAM 
result of 45 zfc 6 K km s _1 . 
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and IRAM data of lGolla fc Wielebinskil l|l994bl ) at the same 
resolution also shows good agreement in the region of the 
central molecular ring. Both the BIMA map and our JCMT 
CO (3-2) map at 17" resolution were then cut off at a con- 
servative 5(7 level before forming the ratio. The resulting 
CO(J=3-2)/CO(J=l-0) map, i? 3 _ 2/1 _ , shown in Fig. [St, 
could be formed only over the brightest, inner 1.8' diameter 
(4.7 kpc) part of the central molecular ring, i.e. approxi- 
mately over the FWHM of the emission shown in Fig. [3] 
The matching smoothed JCMT CO(J=3-2) map is shown 
in Fig. [6]d. 

To estimate the uncertainty in the i? 3 _ 2 /i_o map, we 
have formed a relative error map (not shown), based on 
random errors in the two individual maps. The average of 
the absolute values of the relative errors, is 11% (lower in 
regions of higher S/N). Errors due to positional offsets are 
much lower than this. Allowing for an additional absolute 
calibration error due to uncertainties in t)mb (Sect. |5J, we 
estimate that the average uncertainty in Fig. [6^ is of order 
25%. However, the calibration error can be neglected when 
examining the distribution of _R 3 _ 2 /i_o over the map since it 
would shift all points the same way. We hav e verified that our 
value of -R 3 _ 2 /i-o agrees with the result of llsraell (|2009h for 
the central position when smoothed to the same resolution. 

The average value of i? 3 _ 2 /i_o over the region displayed 
in Fig. [6^ is 0.47 with an rms of 0.11 and extrema of 0.24 and 
0.92. Ninety-five percent of all pixels lie within the range, 
0.28 to 0.66. These variations in _R 3 _ 2 /i_o appear to be real 
since they exceed the approximately 11% uncertainty dis- 
cussed above. 

Lower values of i? 3 _ 2 /i_o, on average, are generally 
found for galaxy disks or for galaxies globally, for exam- 
ple those that are typically found in t he molecular m edium 
along the Milky Way disk (« 0.4, ISandersI 1 19931 ). the 
disks of galaxies, or samples of nearby galaxies (0.2 to 0.7, 
lMauersbergerlll999l : I Wilson et aD 20091. Higher values, on 
the other hand, are seen in the central regions of galax- 
ies (e.g. 0.7 to 1.2 for M 83, 0.9 for M 82, 0.6 for M 51 
iMuraoka et~aHl2007l ; iTilanus et ai1ll99ll ; llsrael et alj|2006t 
respectively), in high excitation regions or regions associated 
with star formation (e.g. up to 1.6 f or the Antennae and 
up to 1.2 for spec ific regions in M 33, iPetitpas et"ai1l2005l ; 
iTosaki et alj|2007l . respectively), or in galaxies in which tem- 
perat ures and/or molecula r gas densities are higher, on av- 
erage ( Maucrsbcrgcr 1999). Therefore, the line ratios found 
in the central molecular ring of NGC 4631 appear to be typ- 
ical of galaxy disks rather than of regions containing strong 
starbursts. 

For the centr al 21 arcsec diameter region of NGC 4631 
only, llsraell \200$ ) finds that the molecular gas is best rep- 
resented by two molecular components, with most of the 
mass (80%) in a cold (T kin — 10 K), tenuous (tiH 2 « 300 
cm~ 3 ) component, and a lesser amount (20%) in a warm 
(Tkin = 150 K), moderately dense (n_H 2 as 10 3 cm -3 ) state. 
Fig. [6^l has now extended the measurement of _ R 3 _ 2 / 1 _ n ove r 
a much larger region than that measured by Israel (2009). 
With a single line ratio, we cannot place strong constraints 
on the state of the molecular gas over the extended region. 
In general, however, if a single physical state is present, 
the relationship between kinetic temperature, molecular 
gas density, and i? 3 _ 2 / 1 _ can be represented by Fig. [7] 
in which we show a large velocity gradient (LVG) model 
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Figure 6. Comparison of emission in the brightest part of the 
central molecular ring. All maps have been smoothed to 17" res- 
olution, cut off at the 5rr level, and are represented in both con- 
tours and greyscale with darker representing higher values. The 
star marks the galaxy's center, (a) The integrated intensity ra- 
tio map, -R 3 _ 2 /i_q. Contours are at 0.35, 0.45, 0.55, 0.65, and 
0.75. (b) The JCMT CO(J=3-2) map in greyscale and contours 
of the inner region of Fig. [T] Contours are at 11, 22, and 33 K 
km s~ 1 and the peak is 48.9 K km s —1 . (c) Molecular hydrogen 
surface mass density, S// 2 , formed from the CO(J=1-0) total in- 
tensity map, adopting X = 2.0 X 10 2t) mol cm -2 (K km s -1 ) -1 . 
Contours are at 125, 200, and 350 Mq pc — 2 and the peak is 380 
Mq pc — 2 . (d) Radio continuum emission from the FIRST sur- 
vey. Contours are at 3, 6, 9, 12 and 18 mjy beam" 1 and the peak 
is at 21.4 mjy beam -1 . The transformation to T,sfr is given in 
Eqn.rj (e) SFE, with contours at 4.1, 5.5, 8.4, and 12 X 10" 10 
yr -1 . The peak is 17.3 X 10~ 10 yr" 1 . 
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Figure 7. Plot of the ratio, R- 



3-2/1-0 



as a function of molec- 



ular hydrogen density, nu 2 (cm -3 ) for a range of temperature 
(labelled) from an LVG code, assuming an abundance per unit 
velocity gradient of Xi 2 qq / (dv/dR) = 10~ 6 pc (km s~ 1 ) — 1 and 
a single component model. 



(e.g. iGoldreich fc Kwanll 19741 ; llrwin fc Aven1ll992l ) with an 
abundance per unit velocity gradient of Xi2 CO /(dv/dR) = 
10~ 6 pc (km s" 1 )" 1 (the latter value as suggested by the 
results of lzhu et af]|2003h . Values of i? 3 _ 2 /i-o from 0.28 to 



0.66 would then imply gas densities of order 10 cm" over 
the range of temperatures showrF^l. Note that this density 
refers to the individual cloud densities that are responsible 
for CO excitation, rather than mean densities within the 
beam. If two components are present throughout the region 
wit h a lower de nsity component dominating, such as found 
by llsraell (|2009T ) for the center, then 10 3 cm" 3 should be 
an upper limit to nn 2 - Therefore, again we find that the 
conditions within molecular clouds in the central region of 
NGC 4631 appear to be typical of low density molecular gas 
regions in galaxy disks (i.e. nu 2 < 10 4 cm" 3 ) rather than 



the <; 10 4 cm 3 which are more t 


ypical of central starburst 


regions (see, e.g. 


Weifi et alJ 12001 


; Hailev-Dunsheath et al.l 


|200S); llono et al. 1 


2007). 



3.4 Molecular Mass and Gas/Dust Ratio 

The distribution of molecular hydrogen mass surface den- 
sity, Eh 2 , in the central molecular ring of NGC 4 631 is 
shown in Fig. \Ejp. This map was formed from the iRandl 
(20003) CO(J=1-0) map and therefore shows essentially the 
CO(J=1-0) distribution. We have used a standard and con- 
stant CO(J=1-0) integrated intensity to H2 conversion fac- 
tor of X = 2.0 x 10 20 m ol cm" 2 (K km s" 1 )" 1 , consistent 
with I Wilson et all (|2009l ). corresponding to 3.2 M pc" 2 (K 
km s" 1 )" 1 of molecular hydrogen (uncorrected for heavier 
elements). The mass distribution shows similar peaks and 
minima as the CO(J=3-2) map; differences are highlighted 
by the ratio map Fig. HJi. For this value of X the total 



15 If Xi2nn I (d v I dR) changes by an order of magnitude (see 
IZhu et alJ l2003). then over the range of interest, the resulting 
density changes by less than a factor of two. 



molecular hydrogen mass over the region shown in Fig. (6ji 
is Mh 2 = 7.4 x 10 8 M Q . 

For this galaxy, there have been several measure ments 
of the value of X from independent line ratio analyses, llsraell 
|2009i ) finds X = 0.3 x 10 20 mol cm" 2 (K km s" 1 )" 1 within 
the central 21 arcsec and iPaglione et all |200ll ) finds X = 
0.5±l ° x 10 20 mol cm" 2 (K km s" 1 )" 1 within the central 46 
arcsec and X = 1.8±i J x 10 20 mol c m" 2 (K km s" 1 ) " 1 out - 
side of the central 46 arcsec. The twojPaglionc et all (|2001T ) 
values agree with each other and also with ou r adop ted value 
within uncertainties. If we adopt the llsraell |2009r i value of 
X for the central 21 arcsec only, then Mh 2 reduces by ap- 
proximately 11% which is not significantly different from the 
above result, given the uncertainties. Although adopting a 
lower central value of X does not significantly perturb our 
calculation of total mass obtained from Fig.[6j;, there would 
be changes in the appearance of this map, should X vary 
with position. 

In addition, since we have observed CO(J=3-2) over a 
larger region than shown in Fig. [6] and to larger radii than 
previously detected in any CO transition, we can estimate 
the total molecular hydrogen mass in NGC 4631 by applying 
an appropriate value of i?3_2/i_o over the entire emission 
shown in Fig. [Tj and using the standard value of X listed 
above. Adopting the mean value from the central region, 
i? 3 _2/i_o = 0.47, we find a total mass of Mh 2 = 2.2 x 10 9 
Mq with an error of « 25% which represents the uncertainty 
in Fig. [6^l including the calibration error. This uncertainty 
does not include uncertainties in X or in its possible varia- 
tion with position in the galaxy. Our result has improved on 
thatofMflv, = 1.5 x 10 9 M fl (adjusted to our dista nce and 
value of X) provided by iGolla fc Hummel (1 199481 ). whose 
map does not extend as far out as ours and who suggested 
a factor of 2 uncertainty on their quantity. 

The total molecular gas mass is 22% of the to tal HI mass 
of M H i = 1.0 x 10 10 M found bv lRandl l| 19941 ) (Table [J). 
The total HI + H 2 mass is therefore M H i+h 2 = 1.22 x 10 10 
Mq and dominated by HI. Adjusted for heavy elements (a 
factor of 1.36), the total gas mass is then M 9 = 1.66 x 
10 10 M The total dust m ass in NGC 4631 i s estimated to 
be M d = 9.7 x 10 7 M Q (|Bendo et al.ll2006l ). leading to a 
global gas-to-dust ratio of 170, a value t hat is typical of spira l 
galaxies, including the Milky Way (e.g. iDraine et al.l l2007). 
These masses are integrated over the entire galaxy and do 
not necessarily represent the relationships between atomic, 
molecular, and dust components in individual regions; we do 
not have sufficient information to determine region-specific 
quantities without a model for each of those components in 
this edge-on galaxy. 

3.5 Star Formation 

In Table [1] we list two estimates of the global SFR, the 
first from the FIR luminosity, and the second from the 
Ha emission corrected for extinction using A 24 /im data 
(Hocrr). The two values differ by about a factor of two. 
SFRs can be determined from a variety of different trac- 
ers but in a galaxy as edge-on as NGC 4631, optical depth 
effects can become large, uncertain, and can vary in an irreg- 
ular fashion with position. The empirical relation for deter- 
mining Ha corr , for example, although considered relatively 
robust, has not been determined for galaxies that are edge- 
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on |Calzetti et abll2007l ). Fig [5] also confirms that there are 
large differences between the shape of the Ho curve and 
other tracers that are not as badly affected by extinction. 
To probe the spatially resolved SFR, then, Ha corr cannot 
be used with full confidence (see also Footnote 1 1711 and the 
FIR luminosity does not have sufficient resolution 16 !. 

An alternative is to use radio continuum emission which 
requires no correction for optical depth effects and for which 
we have data for the high signal-to-noise region at the center 
of the galaxy as shown in Fig. [SJl, taken from the Faint Im- 
ages oFth^Fiadjo Skyat Twenty-Centimeters (FIRST) sur- 
vey (|Becker et al,|[l995h and smoothed from an original spa- 
tial resolution of 5 arcsec. The FIRST survey is insensitive to 
structure greater than 2' in scale which is approximately the 
diameter of the region shown in Fig. HJ] and therefore there 
should be no missing flux on the scales that we are prob- 
ing. In addition, although CR electrons diffuse from their 
source of origin in comparison to other tr acers of SF, the 
scale over which this occurs for NGC 4631 (|Marsh fc Heloul 
1 1998T ) is less than the beam size of Fig. [6] Finally, there is 
no evidence for a radio emitting AGN (or candidate) in this 
galaxy at a flux level that could affect the SFR determina- 
tion (see lGollalll999l ). The radio continuum is therefore a 
good measure of massive SF in the displayed region from 
which we can then estimate the total SFH 17 l 

We first use the relation given in ICondonl i|l992l ) which 
provides only the massive (M > 5 Mq) SFR, SFR m from 
the radio emission, i.e. 



SFR„ 



[M yr- 



= 2.47 x 10" 



Li, 



[W Hz- 



(1) 



where Lia is the 1.4 GHz radio luminosity per unit band- 
width and we have assumed that the 1.4 GHz emission is 
dominated by the non-thermal componeni 18 ! with a typical 
spectral index of 0.8. 

Integrated over the region displayed in Fig. [6jl, we find 
SFR m = 0.34 M Q yr _1 (or v = 0.014 SNe yr _1 , from re- 
lations in Condon, 1992). We can also obtain SFR m for the 
entire galaxy (not just the region shown in Fig. [6]) from 
Eqn. [T] and the globa l radio continuum flux (771.7 mjy 
IStrickland et al.ll2004bf ) resulting in SFR m = 1.8 M Q yr _1 . 
Therefore, « 19% of the massive SF is occurring in the re- 
gion shown in Fig. Thus, the massive star formation in 
this galaxy is widely distributed in contrast to 'nuclear star- 
bursts ' such as M 82 (see also arguments in lTiillmann et all 
l2006bl ). The Galaxy Evolution Explorer (GALEX) UV im- 
age also reveals widely distributed UV emissi on consistent 
with distributed star formation in this galaxy (|de Paz et al.l 
I2007D . 

We can now apply a factor to account for non-massive 
star formation, i.e. SFR/SFR m = 2.4, where SFR is taken 
to be SFRf/j? from Table Q] and SFR m is given above. Note 



16 The FIR data have spatial resolutions of 1.44 a rcmin and 2.94 
arcmin at A 60 /^m and A 100 (mi, respectively jSanders et al.l 
l2003f) . 

17 A map of Ha corr should resemble Fig. [5)1 if the former has 
been adequately corrected for extinction. We have verified that 
there are sufficient differences between the two maps and therefore 
Fig.[6ji is the p referred tracer of spatially resolved SF. 

18 iGollal jl999l ") estimates a thermal fraction of 10% at 5 GHz, 
indicating that the fraction will be lower at 1.4 GHz. 



that this correction factor is mid-range between the correc- 
tio n factors that w ould result by applying a disk IMF given 
bv lChabrierl (|2003T l (a factor of 1.1) and a Salpeter IMF (a 
factor of 5.5) over a total mass range of 0.1 ^ M ^ 100 Mq 
yr -1 . The result, which converts the values of Fig. [6ji to 
SFR per unit area, is 



[M yr" 1 pc- 



= 9.2 x 10 



-6 h.4 

[Jy beam - ] 



(2) 



We now provide a measure of the efficiency 
of star formation, SFE, in the central regions of 
NGC 4631 by forming the ratio map of SFE/[yr -1 ] = 
Ssfr/[Mq yr -1 pc _2 ]/£h 2 /[M0 pc -2 ], which assumes that 
the stars mainly form in molecular gas (excluding HI). The 
result is shown in Fig. Note that we make no correction 
for inclination, so the values in this figure represent SFEs 
over the line of sight through the region. The mean value is 
SFE = 6.4 x 10 -10 yr -1 with an uncertainty of 2.5 x 10~ 10 
yr -1 (based on the rms of the map) and extrema of 3.5 and 
17.3 x 10 -10 yr " 1 . This result fall s withi n the la error bar 
of the sample of iRownd fc Yound (1 19991 ) (after conversion 
to their definition of SFEl 19 !) who examined global SFEs for 
568 galaxies; that is, although there is a strong concentra- 
tion of molecular gas in the central region of NGC 4631 and 
the galaxy is interacting (Sect.[l|, the SFE in this region is 
typical of galaxy disks in general. 

The inverse of SFE is a simple measure of the gas 
consumption timesc ale, t a , if the SFR rem ains constant 
with time. Following iKnapen fc James! (2009) and references 
therein, including a correction for recycling of material by 
stars into the ISM yields t g = 1/(0.6 SFE). For the region 
shown in Fig. [6J5, we find a mean timescale of m 2.6 x 10 9 
yr (to within approximately a factor of two, given the above 
uncertainties). This result will be a lower limit to the total 
gas consumption timescale since it does not include HI. From 
the arguments of the previous section, we expect the mean 
value of t g to increase by a factor of 2 or 3, should HI be in- 
cluded. All uncertainties considered, this r esult is still consis- 
tent with values found i n other gala xies (Knape n fc Jamesl 
Bigiel et al.ll2008l; iLerov et ail l200alKennicutt etld] 



2009 
1994 



Golla fc Wielebinskill994bh . Clearly, the central re- 



gion of NGC 4631 contains a strong build-up of molecular 
gas (Fig. [3} , but the gas consumption timescale is long for 
a constant SFR. We return to this point in Sect. [4] 

There is clearly a similarity between the 7?3_2/i_o map 
(Fig. [6K) , representing molecular gas excitation (higher den- 
sity and/or higher temperature as indicated by Fig. [7)l and 
the SFE map (Fig. [6^) representing the star formation rate 
per unit molecular gas mass. Regions of lower ratio, as noted 
in Sect. 13.41 approximately correspond to r egions of lower 
SFE - a trend also observed in M 83 by iMuraoka et ail 
i|2007h . There are still differences, however. For example, a 
map of the ratio of SFE/i? 3 _ 2 /i_o (not shown) results in an 
rms variation of 25%, the most important difference being 
at the nucleus at which the SFE appears to be enhanced in 
comparison to i?3_2/i_o- Since both of these maps have been 
formed by normalization with the CO(J=1-0) distribution, 



9 The transformation given in iKennicut 3 (|l998f ) was used to 
obtain an Ho SFR but no inclination correction is made. 
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Figure 8. Global profile formed from the 3 central frames of 
Fig. [9] The rms is 0.5 K. The velocity is heliocentric. 



the nuclear enhancement can be easily seen by directly com- 
paring the CO(J=3-2) distribution of Fig. [6p (from which 
J?3_2/i_o has been formed) to the 20 cm radio continuum 
map of Fig. [6ji (from which SFE was derived). The radio 
continuum map shows a strong nuclear peak whereas the 
CO (3-2) map does not. Thus, there is an enhancement in 
SFR and SFE right at the nucleus in comparison to the 
surrounding region. 



3.6 The Velocity Distribution 

The global profile and position-velocity (PV) slices along 
and parallel to the major axis of NGC 4631 are shown in 
Fig. [8] and Fig. [9l respectively. 

It is unusual to show a CO(J=3-2) global profile for a 
galaxy and Fig. [8] serves to illustrate the high quality of the 
HARP data. From this plot, the line width at half-maximum 
is W50 = 270 ± 10 km s _1 and at 20% of maximum is 
W20 = 316 ± 10 km s" 1 . From the mean midpoint of the 
linewidths, we find V sys = 607 ± 10 km s~ , in good agree- 
ment wit h the HI val ue of Table Q] as well as the HI value 
found by iRandl |l994|) (610 km s -1 ) t hough it is somewhat 
lower t han the CO(J=1- 0) value of Golla fc Wielebinska 
il994bT ) (628 km s" 1 ) and lGao fc Solomon! (|2004h (658 km 



This global CO(J=3-2) value of V aya is plotted in 
Fig. [TOji. 

The PV velocity distribution shown in Fig. [5] will be 
considered in three parts, namely, the central molecular ring 
and nucleus, the outer disk emission (both discussed below) 
and anomalous features and high latitude gas (to be dis- 
cussed in Sect. 13. 7l) . 



3.6.1 The Central Molecular Ring and Nucleus 

The most dominant emission shown in Fig. [9] is, again, the 
central molecular ring which extends 3.8' (10 kpc) in di- 
ameter (as noted in Sect. 13. 1[) forming a strongly emitting 



All values corrected to our velocity definition, where necessary. 
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Figure 9. PV plots parallel to the major axis averaged over a 
width of 11" (contours on greyscale). Vertical offsets (positive to 
the north) are marked on each frame. Contours are at 0.044 (2a), 
0.07, 0.10, 0.16, 0.25, 0.35, and 0.6 K. Positive is to the east on 
the x axis (50" = 2.2 kpc). F1-F5 are discussed in Sect. 13. 7. ll 
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region with a steeply rising rotation curve. The velocity gra- 
dient is 2.1 km s -1 arcsec -1 , or 48 km s -1 kpc -1 with an 
estimated 25% uncertainty depending on where the slope is 
measured. The fact that the western peak is stronger to the 
south (offset of -11") and the eastern peak is stronger to 
the north (offset of +11") reflects a slight asymmetry in the 
curvature of the major axis. (The major axis position angle 
of 6.5 deg was adopted to pass through both eastern and 
western maxima.) These gradients ag ree with those of the 
CO(J=1-0) distribution (|Randll2000af ). 

An interesting new result is the 'kink' in the rotation 
curve approximately at the infrared nucleus of NGC 4631 
such that in a single pixel (0.16 kpc) there is a vertical drop 
in velocity (east to west) of 25 km s -1 . Bright emission 
around the nucleus has been emphasized in the 'blow-up' 
of this region shown in Fig. 1 10b ,. The kink is seen at 640 km 
s _1 just to the east (left in the figure) of the nucleus and the 
vertical drop continues to the west of the nucleus where the 
emission is much fainter since it falls within the emission gap 
denoted 'C in Fig.[T] Similar behaviour can be seen in Offset 
+ 11" frame of Fig. [5] and thi s steep gradient is also hinte d 
at in the CO(J=1-0) data of lGolla fc Wielebinskil (|l994bl V 
The nuclear gradient is delineated by the green line in the 
figure which has been adopted to pass through the small 
emission peak at the nucleus (this peak was pointed out in 
Sect. l3.L2j) . Our modelled ring center (from Table[T]) and our 
value of V S ys (from the global profile of Fig. [S]) are marked 
by the cross and agree, within uncertainties, with the center 
adopted for the nuclear gradient. This location appears to 
mark the galaxy's true nucleus and is an additional argu- 
ment (see also Sect. l3.L2j) for adopting the IR center as the 
location of the nucleus rather than the minimum at C. 

The slope of the nuclear gradient is 4.1 km s -1 arcsec - 
(94 km s -1 kpc -1 ), much steeper than the rotation curve 
of the central molecular disk, implying the presence of a 
centrally concentrated mass of Md yn — 5 x 10 7 Mq within 
a radius of 282 pc. To our knowledge, this is the first evidence 
for a concentration of mass at the center of NGC 4631. 



3.6.2 The Outer Disk and Total Dynamical Mass 

The weak outer disk emission becomes most evident at radii 
greater than 75" (3.3 kpc) and appears kinematically dis- 
tinct from the central molecular ring (see Fig. [9]). Peak ve- 
locities at the largest measurable radii on either side of the 
nucleus are approximately the same as the peak values found 
in the central molecular ring, but in the region between the 
ends of the central molecular ring and the furthermost radii, 
velocities are lower, giving the impression that the rotation 
curve declines at the ends of the central molecular ring and 
then rises again with radius. However, a comparison with 
PV plots from the HI distribution shows that the rot ation 
curve does not decline in this region (see iRandl |l99J) , his 
Fig. 9); rather, it simply becomes flat at the ends of the 
central molecular ring and outwards. Therefore the appear- 
ance of lower rotational velocities just outside of the central 
molecular ring is a result of irregularities in the CO(J=3-2) 
emission intensity (see next section). Given the faintness of 
the CO(J=3-2) emission in the outer disk, it has not been 
possible to identify features corresponding to the HI super- 
shells in these regions (Fig. [I}. 

Taking the mean of the maximum velocities (V r = 155 



km s -1 ) and radii (R = 10.6 kpc) on either side of the 
nucleus and adopting a spherical distribution of total (light 
plus dark) mass for NGC 4631, the total mass is M to t(R < 
11 kpc) = 6 x 10 10 M©. The HI distribution reaches compa- 
rable velocities (150 km s -1 ) but ca n be detecte d to much 
larger radius, i.e. to R « 24 kpc (|Randl fl99i ) . From HI 
data, we can therefore extend the total mass estimate to 
M tot (R < 24 kpc) = 1.3 x 10 11 M . 



3.7 Anomalous Velocity Features and High 
Latitude Molecular Gas 

3.7.1 Anomalous Velocity Features 

Fig. [9] reveals several anomalous velocity features which ap- 
pear as extensions or partial loops in PV space associated 
with the central molecular ring. We have identified five such 
features (labelled Fl through F5 on the major axis slice) 
labelling only those that are connected to the central molec- 
ular ring and can be traced over at least two beam sizes spa- 
tially, at least two contiguous velocity channels, and over 2a 
in intensity. Two features, Fl and F3, occur near the ends 
of the central molecular ring and contribute to the (inac- 
curate) appearance of a declining rotation curve in these 
regions. To aid in visualizing these features, we have spa- 
tially smoothed the major axis slice of Fig. [5] and show the 
result in Fig. 110b with the same labe lling. The inset shows 
the CO(J=1-0) data of lRandl (^OOOaT ) treated similarly. Al- 
though the CO(J=1-0) data show significant noise, they do 
reveal some extensions at approximately the same locations 
of the features we have identified. We have also verified that 
these features exist in the CO(J =3-2) data cube tha t has 
been independently reduced by IWarren et all l|201Ch (see 
Sect.©. 

Some of the labelled features can be traced above or 
below the major axis. For example, the feature, F3, which 
appears loop-like on the major axis slice of Fig. [9] shows a 
split velocity profile above the major axis (Offset = +11") 
open to the east. This feature c orresponds to the expanding 
shell observed bv lRandl |2000al ) in the CO(J=1-0) distribu- 
tion. The feature, F2, shows complex structure (Tig. I10p but 
it can be traced south of the major axis (Offset = -11" in 
Fig. [9| where it has the appearance of a smaller complete 
loop. The feature, Fl, could be related to F2, though this 
association is not clear. Feature F4 appears to be an anoma- 
lous velocity extension to the main emission, and F5 forms 
a large, but weaker loop. 

Since these features blend with the main emission in 
RA-DEC space, we cannot confirm shell-like structure spa- 
tially; however, the appearance of at least F2 and F3 
in PV space are consistent with known behaviour of ex- 
panding shells or portions of expanding shells such as 
have be en seen in our own and other galaxies in HI or 
CO fe.g.lSpekkens et al.ll2004)McClure-Griffiths et al.ll2002l ; 
llrwin fc Sofudll996l ; lLee et al(l2002l ; IWeifi et al.lll999ly F3 is 
the first detection in CO(J=3-2) of the previously known 
CO(J=1-0) shell and the other features are new detections. 

In Table[4]we provide parameters of these features. The 
molecular masses may be lower limits because CO(J=3-2) 
emission represents only gas that shows some excitation and, 
in addition, we have determined the masses only over regions 
where the features clearly depart from the main disk emis- 
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Figure 10. Enhanced versions of the major axis frame of Fig. [9] (a) Blow-up of the nuclear region emphasing the brightest emission 
with contours at 0.2, 0.25, 0.3, 0.33, 0.35 and 0.4 K. The x axis is centered on the IR center (Table[T}. The cross marks the location of 
the modelled ring center in x (Table [3} and V sys (from the total global profile of Fig. |5J, with the cross size delineating their respective 
error bars. The green line shows the nuclear gradient (see Sect. 13. 6. It . (b) Here, the major axis has been smoothed spatially with a 
Gaussian of FWHM = 11" to emphasize the fainter emission. Contours are at 0.028 (2<r), 0.035, 0.050 and 0.074 K. Features have been 
labelled as in Fig. [9] Inset: CO(J=1-0) major axis emission from the original data of Rand (2000a) rotated and smoothed as in the main 
figure; contours have been arbitrarily set. 



sion in PV space. Nevertheless, we do find that our m ass for 
F3 agr ees with the value of « 10 8 Mq measured by iRandl 
|2000al ) for the corresponding CO(J=1-0) expanding shell. 

The velocity width of these features, A V, indicates 
that, even if no shell-like signature is observed, some ex- 
pansion Vexp = AV/2 must be occurring. From this in- 
formation, we estimate the kinetic energy of the anomalous 
velocity features, Ek (see Table S}. These energies may also 
be lower limits since they scale with mass. Although we are 
not certain of the origin of the features (see Sect. 3), it is 
common practice to estimate the input mechanical energies 
required to form them, Eo, assuming an instantaneous en- 
ergy input, such as would be the case for clustered super- 
novae and stellar winds over a time scale that is short in 
comparison to the age of t he feature itsel f. The required 
mechanical input energy is l|ChevalierJll974l ) . 



= 5.3 x 10 43 
ergs V cm 



P4) ] 
V cm _J / 



Re 



pc 



km s 



(3) 

where no is the ambient density at the time of the energy 
deposition and R exp — D/2 is the shell radius. If a contin- 
uous energy input is assumed instead, the results tend to 
be consistent with E qn. [3] to within a factor of a few (see 
ISpekkens et all2004l ) . Models such as these typically assume 
a uniform ambient density in the plane which is certainly not 
the case for molecular clouds. For NGC 4631, most of the 
molecular mass is in clouds with densities ~ 300 cm -3 (see 



Sect. 13. 3[) . Since we estimate that there are roughly equal 
HI and H2 masses in the regio n of the cen tral molecular 
ring (Sect. 1374)) and the model of lRandl l|l994h estimates the 
ambient HI density in the plane to be 0.4 cm -3 , then the 
average H2 density could be similar, leading to a volume 



filling factor of ?a 10 -3 for the molecular clouds. Although 
these are very rough estimates, most of the molecular out- 
flow would be through in-plane densities that are lower than 
those of the molecular clouds themselves. We have therefore 



used no 



0.4 cm 3 in Eqn. [3j bearing in mind that the 



results for Eq will be conservative estimates since higher 
densities could be important at some level. 

Finally, we compute 'characteristic' ages, r, which do 
not take into account accelerations or decelerations over the 
development of the feature. If a continuous wind model is 
adopted, these lifetimes would decrease by ap proximately a 
factor of three (|McClure-Grifnths et alj|2002 ). It is interest- 
ing that all lifetimes fall into a narrow range of timescales, 
r ~ 2 — 5* 2.6 x 10 7 yr, suggesting that they may be re- 
lated to a single burst of star formation. Note, however, that 
our spatial resolution selects features that are of kpc scale 
and these observations would not have detected smaller, and 
therefore younger, features, if they were present. Larger fea- 
tures may also be difficult to detect if their densities diminish 
with increasing size. 

The CO(J=3-2) shell results of Table H are similar 
to those of HI shells found in our own Milky Way and 
external galaxies (e.g. | Brinks fc Baiaial Il986l: IPuche et alj 



19921; iMcClure-Griffiths et all |2002| ; IChaves fe Irwinl l200ll ; 
Spck kens et al J2004L and others). The tabulated masses and 
energies, although order of magnitude estimates, are likely 
conservative as discussed above; it is clear that many hot 
young stars and supernovae would be required to form the 
features if they are indeed the origin. We will return to this 
issue in Sect. [4] 
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Table 4. Parameters of the Anomalous Velocity Features 



Feature 


RA a , DECo a 
(h m s, ° ' ") 


Vol 

(km s" 1 ) 


c, k^c) 


AV^ 
(km s~) 


M mol ± 
(10 7 M ) 


(10 53 ergs) 


E f_ 

(10 53 ergs) 


t9_ 

(10 7 ~r) 


Fl 


12 42 00, 32 32 15 


500 


22, 0.95 


52 


2.3 


1.6 


4.1 


1.8 


F2 


12 42 03, 32 32 13 


593 


44, 1.9 


73 


5.3 


7.1 


57 


2.5 


F3 


12 42 14, 32 32 40 


718 


33, 1.4 


62 


9.0 


8.6 


18 


2.2 


F4 


12 42 10, 32 32 29 


593 


22, 0.95 


42 


3.0 


1.3 


3.0 


2.2 


F5 


12 42 07, 32 32 25 


468 


44, 1.9 


73 


2.4 


3.2 


57 


2.5 



a Central position and velocity of the feature, from the centroid of the emission after smoothing spatially and in velocity 
using both the major axis slices as well as offset slices (see Fig. [9j as needed. Positional uncertainties are approximately ± 10" 
and the velocity uncertainty is ± 10 km s — 1 . 

b Diameter of the feature (angular and linear) measured from the original resolution data. Uncertainties are as in Footnote [a] 
c Full velocity extent of the feature measured from the original resolution data. Uncertainties are as in Footnote [a] 
^ Total molecular mass (including heavy elements), adopting R 3 _ 2 /i-o = 0.47 (Sect. IXHl ; the result is an average between the 
smoothed and unsmoothed data. The uncertainty is «i ±0.5 X 10 7 Mq which represents a typical flux in the background over a 
similar-sized region. 

e Kinetic energy of the feature from E K = (1/2) M (A V/2) 2 . 
J Input energies, from Eqn.[3] 

3 Characteristic age of the feature, from r = D/(AV). 



3.7.2 High Latitude Molecular Gas 

As indicated in the previous section, the anomalous velocity 
features seen in the PV plots (Fig. [9| are not easily traced 
to high latitudes. However, we do see some evidence for the 
presence of high latitude CO(J=3-2) in the data. 

Fig- H shows that the disk thickness of NGC 4631 varies 
with position, but at least within the central molecular ring, 
the evidence suggests that NGC 4631 forms a thick, rather 
than a thin distribution. A thin global molecular gas disk 
with the observed diameter of the CO(J=3-2) emission (see 
Sect. 13. 672"]) which is inclined by 86° (Table[T]) would project 
to a total apparent vertical extent of only 48" including the 
smoothing effects of the beam, whereas the total observed 
minor axis extent is approximately 65" (2.3 kpc after beam 
correction, or z — 1.2 kpc). If the inclination of the central 
molecular ring were as low as 83°, then it could be inter- 
preted as thin. However, our model of the central molecular 
ring (Table O gives best results for an even higher inclina- 
tion (i — 89°); lower values are poorer. 

To further investigate the vertical extent of the molecu- 
lar gas, we have formed a plot at high sensitivity showing the 
minor axis profile by averaging over a 100"-wide region of the 
major axis (i.e. approximately over the FWHM of the central 
molecular ring) and then averaging the north/south sides. 
The result, shown in Fig. 1111 reveals a complex profile which 
is not well described by a single smooth fit. CO(J=3-2) can 
be traced out to z = 33" (1.4 kpc), correcting for beam 
smoothing. Thus, the central molecular ring of NGC 4631 
appears to have a thick vertical distribution of CO. Given 
the known halo activity in this galaxy, most of which has 
been measured to much larger values of z (see Sect.[T]), the 
presence of thick, agitated molecular gas is perhaps not sur- 
prising. For comparison, we note that main sequence stars 
in NGC 4631 have been measured out to a z height of 2.3 
kpc and AGB a nd RGB stars hav e been measured to even 
higher z values (|Seth et al.ll2005bl ). 

We see no evidence from the PV plots (Fig. [9]), however, 
for a change in the slope with z height ('lagging halos') such 




Distance from Major Axis (arcsec) 



Figure 11. Vertical profile of the CO(J=3-2) emission of 
NGC 4631, formed from Fig. [T] averaged over a 100" wide re- 
gion centered at the nucleus. The profiles to the north and south 
have also been averaged to form a sensitive plot. The emission 
cuts off at z m 40" due to the algorithm for determining the 
total intensity image (see caption of Fig. [TJ. The projected mi- 
nor axis extent of a thin disk of radius, R = 243", inclined at 
i = 86° (Rcos(i)) and then smoothed by the beam, is given by 
the thick vertical line segment. 



as has been seen in HI or ionized extraplanar gas in sev- 
eral other galax i es (e.g. iTiillmann et alll20o"cl ; lRandll2000bl ; 
iFraternali et~ai1 120021 : lOosterloo et al l 120071 ) although the 
vertical extent of molecular gas in NGC 4631 is, in general, 
smaller than these other components. 

Fig.[T]also shows a number of disconnected emission fea- 
tures above and below the plane at distances of 50" — 60" 
(z ~ 2.2 — > 2.7 kpc). The smoothing and noise cut-off tech- 
niques used to create the total intensity image are ideal for 
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emphasizing such low level emission which generally can only 
be seen in cubes that are smoothed, rather than in the orig- 
inal channel maps. Their emission is contiguous in velocity 
space and several can be identified in the independently re - 
duced, lower velocity resolution data of lWarren et al.l |2010). 
However, since not every one can be independently con- 
firmed, we do not label each individually and only consider 
the energetics (below) of a 'typical' feature. For the purpose 
of discussion, we refer to them as high latitude 'clouds'. 

Fig. 1121 which emphasizes the low level emission of 
Fig. [T] shows how these clouds exist above apparent dis- 
turbances in the disk. For example, there are clouds both 
above and below the location of the western HI supershell 
(the western X). There is also a cloud above the plane and a 
complete loop below the plane at the location of feature F5. 
These clouds have masses of M c ; ~ 10 6 Mq and estimated 
potential energies of, 



5M = 2 x 10 " 
ergs 



2 x 10 7 M 



700 pc 



(h) 



x In 



cosh 



(4) 

,O.185M pc~ 3 j [ \z J\ 

(e.g. lSpekkens et al.ll2004T ). where zq is the mass scale height 
and po is the mid-plane mass density with the normalization 
factor taken to be an estimate for the local Solar neighbour- 
hood. Adopting th e stellar scale height of Zo = 1 kpc from 
ISeth et all (|2005bh we find, 

e|s - 3 ' 7 X 1053 ( o.l85M B pc-» ) (5) 

which suggests that kinetic energies of this order are required 
to place the molecular clouds at their observed heights. In- 
put energies are likely higher still (cf. Ek, Eo in Table [3} if 
these clouds have originated in the disk. Further discussion 
can be found in Sect. [4] 

3.7.3 Comparison of Outflow and High Latitude Features 
with other Wavebands 

As indicated in Sect. [T] the halo of NGC 4631 has been ob- 
served in every ISM component. Since the halo is so exten- 
sive and there is overlap with the broader HI tidal streamers, 
some caution must also be exercised in identifying correla- 
tions as they may be due to chance projec tions along the 
line of sight (see e.g. iTavlor fc Wand [2003). The features 
that we have identified (Table 0]) are mostly visible in PV 
space, rather than in RA-DEC and we have thus not found 
counterparts at other wavebands, other than F3 as noted in 
S ect . 1 3 . 7TT1 which is associated with a CO(J=1-0) expanding 
shell. 

A possible hig h latit ude dust arch, discovered by 
iNeininger fc Dumkel (1 19991). bu t inter preted as part of HI 
tidal spur 4 bv lTavlor fc Wand l|2003l ). has its footprint in 
the disk at approximately the position of the west HI ex- 
panding shell. At this location, early X-ray images show 
hot gas extending into the halo, as shown most clearly in 
IWang et ail (|l99a ) (their Fig. la). More recent X-ray im- 
ages (see Fig. I12[) show emission from hot gas which has 
originated from the disk and i s associated with H a emission 
and a 'froth' of superbubbles ([Wang et al.ll200lh . Given the 
varying distributions of the X-ray and CO(J=3-2) compo- 
nents and the edge-on nature of the galaxy, Fig. [12] does 
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Figure 12. CO(J=3-2) map of Fig.[T]with low level emission em- 
phasized overlaid with X-ray contours representing hot gas (from 
the data of Wang et al. 2001). Approximate locations of the fea- 
tures identified in Fig. [9] are marked as well as the two HI shells 
(marked with X) identified by Rand & van der Hulst (1993). 



not show clear correlations between the two components. 
However, it is clear that the stronger X-ray halo emission is 
located over the region of the central molecular ring. 

Finally, one high latitude CO(J=3-2) feature, namely 
the southern loop associated with F5 (Fig. 1121) . does appear 
to have an Ha counterpart in the form of two Ha spurs, the 
latter seen in Fig. [4^. To the east of the two spurs is a third 
Ha spur that appears to have a counterpart in the A 160 fim 
emission (Fig. 2b) • 



4 DISCUSSION 

A picture has emerged of NGC 4631 as an actively star form- 
ing galaxy, but not one with a central starburst. Approxi- 
mately 80% of the massive star formation in this galaxy 
occurs outside of the central 1.8' (4.7 kpc) diameter region 
and Ha and UV emission in the disk are also widely dis- 
tributed. The X-r ay halo (Fig.[T2"l) . which resembles the well- 
known radio halo l|Wang et al.ll200ll ). is similarly widespread 
over the entire star forming disk. Widespread halo emission 
is seen in all ISM components (see Sect. [T] for references) 
and all evidence, including the presence of HI supershells, a 
CO(J=1-0) expanding shell, vertical filaments, the vertical 
structure of the magnetic field lines, the concentration of 
stronger halo emission over stronger star forming regions in 
the disk, and now CO(J=3-2) anomalous velocity features 
point to star forming regions and their related supernovae 
and stellar winds in the disk as the origin of the halo. The 
question is, why does NGC 4631, whose halo is arguably the 
most spectacular ever observed, have such a wide-spread, 
prominent halo in comparison to other edge-on galaxies? 
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To consider this question, it is worth comparing 
NGC 4631 to two other galaxies: NGC 5775, which has 
a similarly w ide-spread multi-phase halo (|Lee et al.l l200ll ; 
iLi et al.ll2008l ). and M 82, which is the prototypical nuclear 
starburst with bipolar outflow. All three galaxies have sig- 
nificant companions with which they are interacting, the 
most obvious evidence be ing the presence of HI bridges or 
streamers between the m l|Yun, Ho. fc Loi ri994; Ilrwinlll994l ; 
IWeliachew et al.lll978T l. For NGC 4631, the interaction is 
with the large spiral, NGC 4656 and the dwarf elliptical, 
NGC 4627 (Sect.[JJ. These interactions may have played an 
important role in forming their halos via tidal disruptions 
or agitation of the disk. For example, if the disk has become 
'puffed up', such as we seem to see in the thick CO(J=3- 
2) layer in NGC 4631 (Sect. 1377^1 . it would be easier for 
material to escape into the halo. Stellar winds and super- 
novae should also produce thickened gaseous disks as ma- 
terial escapes vertically into the halo. For NGC 4631, since 
the stellar disk also extends several kpc above the plane 
(Sects. [1] I3.7.2[) . the interaction has most likely played an 
important direct role. 

As for star formation, M 82, NGC 5775 and NGC 4631 
all have similar FIR luminosities to within a factor of 
about 2, but their luminosities per unit optical disk area, a 
distance-independent qu antity, vary in the ratio , 1:0.16:0.05 
in the order listed above (|Tullmann et al.ll2006al ). NGC 5775 
and NGC 4631, which both show disk-wide halos, differ by 
only a factor of 3. M 82, on the other hand, has a SFR per 
unit area that is higher than that of NGC 4631 by a factor 
of 20. Since M 82 has been well-studied, it is possible to 
compare activity within the nuclear region itself. For M 82, 
the supernova rate in the central 700 pc is v = 0.1 yr - 
jKronberg et al. T1981J ). From Eqn. [TJ Fig. [6jl, and relations 



Condonl jl992i r we find v — 0.001 yr 1 for an equiva- 



lently sized region at the center of NGC 4631. That is, the 
central supernova rate in NGC 4631 is two orders of mag- 
nitude smaller than in M 82, hence the nuclear outflow in 
M 82 but not in NGC 4631. 

It is well known that interactions can trigger a build- 
up of molecular gas in galaxy centers an d can also induce 
strong central starbursts (e.g. lAaltoll2007i ). What is partic- 
ularly striking about NGC 4631 is the presence of strong 
central molecular emission in a ring out to a radius of ~ 5 
kpc (e.g. Fig. [3]) but without a central starburst. The phys- 
ical properties of the ring are typical of galaxy disks, rather 
than other known nuclear starbursts (Sect. [33} an d the star 
formation efficiency is not very high; that is, if SF contin- 
ues at the current rate, there is sufficient molecular gas to 
sustain it for at least 3 x 10 9 yr. Inside the ring, within 
the central 17" (740 pc) diameter region (Fig. [6j the SFR 
is enhanced and there is a peak in the hot dust distribution 
(Fig. [SI as well as a small peak in CO(J=3-2) (Fig. EJ). Even 
here, however, the gas consumption time scale is at least 
10 9 yr. If the interaction is to trigger a central starburst 
in NGC 4631, then it has already happened or it has not 
hap pened yet. 

iKnapen fc James! (|2009h have added further confirma- 
tion to the notion that, although bursts of SF can occur as 
a result of interactions and there is a tendency for the star 
formation to be centrally concentrated, interactions can also 
initiate continuous star forma tion over longer timescales, 
ie. a few x 10 8 yrs (see also |Pi Matteo et al.l [2008). For 



NGC 4631, the inte raction with its companion occurred 
~ 3 x 10 years ago (|Combeslll978T ). There is some evidence 
that at least one burst of star formation (or higher SFR) 
has already occurr ed in NGC 4631. From UV observations 
I Smith et al.l (|200ll ) note that the current rate of SF does 
not seem to account for other indications of strong outflow 
in this galaxy. For example, the eastern HI shell (see Fig. [TJ 
can be explained by supernovae from a massive SF region 
containing 5.3 x 10 4 OB stars beginning about 2 x 10 7 yr 
ago and that the currently observed UV emission (which is 
insufficient to explain the shell) is due to second generation 
stars. Our own estimates for the lifetimes of the observed 
CO(J=3-2) anomalous velocity features are approximately 
the same (Table [4j and , estimated the same way, th e HI 
supershells observed by iRand fc van der Hulstl (| 19931 ) also 
result in ages of a few 10 7 yr. These results suggest that a 
higher rate of SF may have occurred of order w 10 7 yr ago. 
However, as noted in Sect. 13.7.11 selection effects may have 
prevented us from detecting features with larger or smaller 
lifetimes, so we cannot place a limit on the timescale for past 
SF in general. 

Estimates for the kinetic energies of the anomalous ve- 



locity features are of order 



10 ergs and estimates 



of potential energy of the small clouds above and below the 
plane are about the same (Sect. [3^7. 2\ . The implied input 
energies are higher still, possibly by an order of magnitude 
(Table [4]), a conclusion also implied for the energetics of 
the outflowing molecula r gas in the nuclear outflow of M 82 
(|Seaauist fc Clarkll200ll ). These energies are also typical of 
what has been previously seen for HI and CO expanding 
shells in other galaxies (Sect. I3.7.1[) and are quite high. If 
the observed features are associated with a higher SFR in 
the past as suggested above, however, the energies may be 
sufficient. Supernovae from the stars that are believed to be 
responsible for the eastern HI supershell could supply (at 
10 51 ergs each) 5.3 x 10 55 ergs for this shell which is more 
than adequate. However, there would need to be at least 
five such SF regions throughout the region of the central 
molecular ring to account for all of the features of Table [4] 
There is still some need for time-dependent modelling of 
such outflows in a realistic multi-phase, multi-density, and 
magnetised ISM. 

The prominent radio halo in NGC 4631 actually reflects 
an integration over past SF activity in the disk. The av- 
erage magneti c field strength in the disk of N GC 4631 is 
B = 6.5 fiG l|Dahlem. Lisenfeld fc Golialll995l ) with a CR 
electron lifetime estimate of ten = 4.8 x 10 7 yrs. The ra- 
tio between the halo and disk magnetic field strength is 5/8 
(|Hummel et al.1 1991) yielding an average halo field strength 
of B = 4 fiG and, since ten tx B~ 3 ^ 2 , the CR lifetime in 
the halo is ten ~ 10 8 yrs. Consequently, the radio halo that 
we see today has a memory of outflows that have occurred 
since approximately the time that the interaction occurred, 
assuming that the outflows have not escaped into the inter- 
galactic medium. Whether or not this assumption is correct 
requires deeper and more extensive mapping of the mag- 
netic field direction t han is currently available . However, ar- 
guments presented in lSeaquist fc Clark! I|200ll ) suggest that 
even in the more energetic nuclear outflow of M 82, at least 
the molecular gas and dust do not escape. 

The implication of a higher SFR in the past suggests 
that, when strong halos such as in NGC 4631 and NGC 5775 
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are observed, they may have been enhanced by a previous 
outflow event (or events). Although rather speculative, an- 
other possibility is that such galaxies have also experienced 
a stronger M 82-like starburst and nuclear outflow in the 
past that was triggered by the interaction. In the case of 
NGC 4631, the enhanced SFR right at the nucleus could 
be a remnant of a past 'M 82-like' nuclear starburst. We 
could imagine the interior of the nuclear molecular ring to 
have been excavated by massive star formation and outflow. 
The minimum near the nucleus at 'C (Fig. [T]) could reflect 
local variations in SFR and molecular content. The past ad- 
ditional 'boost' of outflow material into the halo may be 
what is required to distinguish between spectacular halos 
and those that are more modest. 



5 CONCLUSIONS 

As part of the JCMT Nearby Galaxies Legacy Survey, we 
have mapped the CO(J=3-2) emission from the edge-on 
galaxy, NGC 4631, which is known for its spectacular multi- 
phase gaseous halo. 

Most of the CO(J=3-2) emission is concentrated within 
a radius of ~ 5 kpc. Although the spatial distribution could 
be more complex, the emission is well modelled by a sim- 
ple edge-on ring with a Gaussian density distribution which 
peaks at a radius of 1.8 kpc and has inner and outer scale 
lengths of 0.1 and 0.28 kpc, respectively. The center of the 
ring agrees with the IR center. A small CO(J=3-2) peak oc- 
curs within this ring, right at the nucleus. Outside of the 
central molecular ring, weaker more extensive disk emission 
is present which has been mapped out as far as 9.25 kpc 
to the east of the nucleus and 12.4 kpc to west. This radial 
extent exceeds that of any previous CO observations. 

Comparisons have been made between CO(J=3-2), 
A 24 fim, A 160 /im emission and Ha. We find that the 
CO(J=3-2) emission more closely follows A 24 /xm (a hot 
dust tracer) rather than A 160 /J,m emission (a cold dust 
tracer), suggesting that CO(J=3-2) is a good tracer of star 
formation. H a emission is uncorrelated because of the high 
extinction in this edge-on galaxy. 

For the inner 2.4 kpc radius region of the central molec- 
ular ring, we have formed the first spatially resolved maps 
of i?3_2/i_o, the H2 mass surface density, and the SFE for 
NGC 4631. Only 20% of the global SF occurs in this re- 
gion. We find that i? 3 _ 2 /i_o is typical of galaxy disks, in 
general, rather than of regions associated with central star- 
bursts. Molecular cloud densities (« 10 3 cm -3 ) in this re- 
gion are also typical of molecular clouds in galaxy disks 
rather than central starbursts. The SFE in this region is, 
on average, 6.4 x 10~ 10 yr -1 leading to a mean gas con- 
sumption timescale of 2.6 x 10 9 yr for the H2 and longer 
if HI is included. That is, the SFR in this region is modest 
when compared to the abundant gas that is present. There 
is, however, an enhanced SFR and SFE right at the nucleus 
(within a central region of 740 pc diameter), although the 
gas consumption timescale is still long (10 9 yr). 

The total molecular gas mass in NGC 4631 is Mg 2 = 
2.2 x 10 9 Mq, which improves upon previous values. Since 
the total HI mass is Mhi = 1-0 x 10 10 Mq, the total gaseous 
content of the galaxy is dominated by HI. The global gas- 
to-dust ratio is 170. 



The velocity field of NGC 4631 is dominated by the 
steeply rising rotation curve of the central molecular ring 
followed by the flatter outer disk; the peak rotational veloc- 
ity is V = 155 km s - . The total dynamical mass within 11 
kpc radius is 6 x 10 10 Mq . At the center of the galaxy, we 
find a steep rotation curve, providing the first evidence for 
a central concentration of mass, i.e. Md yn — 5 x 10 7 Mq 
within a radius of 282 pc. 

We can now add CO(J=3-2) emission to the long list 
of evidence for outflowing gas in NGC 4631. Five anoma- 
lous velocity features with properties similar to those found 
in expanding shells (or parts thereof) in other galaxies have 
been detected, all associated with the central molecular ring. 
One of these (F3) corres ponds to an e xpanding CO(J=1-0) 
shell previously found by Rand (2000a). The galaxy also has 
a thick CO(3-2) disk which we trace toaz height of 1.4 kpc. 
Some small 'clouds' are observed at higher latitudes, possi- 
bly associated with outflows from the disk. We suggest a 
scenario in which interactions with the companion galaxies 
in the past has produced enhanced star formation through- 
out the disk and speculate that there could have been a 
massive nuclear outflow in the past. 
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